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Abstract 

We propose a simple and accurate method to evaluate the outage probability at the output of arbitrarily fading 



' L-branch diversity combining receiver The method is based on the saddlepoint approximation, which only requires 

the knowledge of the moment generating functions of the signal-to-noise ratio at the output of each diversity 
branch. In addition, we show that the obtained results reduce to closed-form expressions in many particular cases 
' of practical interest. Numerical results illustrate a very high accuracy of the proposed method for practical outage 

'NT \ 

£^ ■ values and for a large mixture of fading and system parameters. 
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Cooperative Communications, Diversity Combining, HARQ, Hoyt fading, Nakagami-m fading. 
Outage Probability, Rice Fading, Saddlepoint Approximation. 

I. Introduction 

IN THIS work, we propose a simple method to evaluate the outage probability at the receiver after 
combining independently fading signals whose distributions are known. Outage analysis is a funda- 
mental problem in communications theory, and has been extensively studied in the literature. It consists 
in finding the cumulative distribution function (CDF) of the sum of independent random variables. Exact 
expressions (e.g., [[D-[|4l) or closed-form approximations Q can be obtained in particular cases and, in 
more general situations, the problem can be solved using a numerical integration (via the inverse Laplace 
transform) [6|. 

The work was supported by the government of Quebec, under grant #PSR-SIIRI-435. When this work was submitted for publication 
M. Jabi was at INRS, Montreal, Canada. 



2 

The approach we propose is valid when fading distributions are arbitrary but with known moment 
generating functions (MGFs). We propose to use the so-called saddlepoint approximation (SPA) [71, a 
well known tool of statistical analysis that was already applied to solve various problems in the area of 
communications [[8l- lfT3ll . 

Using SPA, we are able to approximate — very accurately — the outage probability of diversity combining 
receivers over arbitrarily fading channels. The problem at this level of generality was solved in [l6l using 
Laplace transform. Here, unlike [6], we are able to provide closed-form solutions in two particularly 
interesting cases: (i) when combining L identically distributed signals (corresponding to Nakagami-m, 
Rice, or Hoyt fading), and (ii) when combining L = 2 non-identically distributed Nakagami-m fading 
signals. Our closed-form approximations of the outage probability are accurate and very simple, while 
calculating the exact expressions requires numerical integration, as only in particular cases, it is known 
in closed-form. 

IL System Model 

We consider a communication system where the same signal is received on L independent diversity 
branches at the receiver. While the term branches evokes the receiver connected to different front-ends 
(which reflects the combining implementation in multi-antenna receivers), the scenario we consider also 
covers relay-based communications and automatic repeat request (ARQ) retransmissions. In these cases 
the l-th diversity branch will thus correspond, respectively, to the l-th relay's transmission and the l-th 
ARQ transmission round. 

We assume that the channel between the transmitter and the receiver is varying (fading) randomly from 
one transmission to another but stays invariant during each of the transmissions (block-fading channel), 
thus the signal received on the l-th diversity branch is given by 

y« = a/tI-x + zz, (1) 

where z; is a zero-mean, unitary variance Gaussian signal modelling noise, x is the unitary-variance trans- 
mitted signal, and 7; is the instantaneous signal-to-noise ratio (SNR) on the l-th branch. We assume that 
7; is a random variable whose probability density function (PDF) is given by p^^ (7) and the corresponding 
MGF - by Aii{s) = E^Je'^' '^], where Ex[-] denotes the mathematical expectation calculated with respect 
to X. The average SNR is denoted by 7^ = E^j7i]. In this work, we consider three of the most common 
fading distributions, namely, Nakagami-m, Rice, and Hoyt. Table [I shows the corresponding MGFs. 
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Table I 

MGF Mi{s) OF THE SNR and the derivatives of the corresponding CGF Ki(s) = log Mi{s) FOR THE ADOPTED FADING 

MODELS. 



At the receiver, maximum ratio combining (MRC) of the L branches is performed, and the combined 
signal can be written as 

L 

y = X] V^iyi = ■ ( V7 ■ X + z) , (2) 



n=l 



where 7 = J2t.=i 7' '^'^^ aggregate SNR after combining and z is a zero-mean unitary variance Gaussian 
signal modelling the equivalent noise at the output of the combiner. 

Let 7th be the SNR threshold at the receiver below which the communication is in outagqj. The outage 
probability can thus be expressed as 

Pout = Pr{7< 7th} =F-, (7th), (3) 

where F^(x) = /Q^p^(7)d7 is the CDF of 7. Clearly, finding the outage probability boils down to 
calculating the CDF F^(x). 

Using an infinite- series representation [[T], the CDF can be obtained for Nakagami-m distribution for 
arbitrary values of 7^ and shape parameters uii, cf. Table HI Also for Nakagami-m with integer parameters 
nil, closed-form expressions exist in selective diversity combining scenarios flU. 

An altemative approach to evaluate (|3]) is based on the inverse Laplace transform [|6l. Since 7 is the 
sum of L independent random variables, the Laplace transform of its PDF (or, in other words, its MGF) 
is the product of L individual Laplace transforms. As explained in [|6l, using the MGFs of the variables, 
the outage probability in dD can be found using numerical integration. 

Here, we propose an approach that keeps the generality of the solution in [O valid for arbitrarily 
distributed variables 7^, as long as their MGFs are known. However, our solution (i) does not require 

'Xiiis could be seen as the minimum SNR equivalent to an error rate corresponding to a QoS constraint or, in the case of coded transmissions, 
to the SNR below which successful decoding is not possible 



numerical integration used in [|6l, and (ii) yields closed-form solution for identically distributed 7;. 
Moreover, using our approach, a closed-form outage expression is obtained in the case where L = 2 
Nakagami-m distributed signals are combined. The solution is valid for any values of nii and ~i, and it 
is an altemative to the analytical results shown in 01 which are limited to integer mi. 

III. Saddlepoint Approximation of Outage Probability 

We use here the so-called saddlepoint approximation which is a simple and accurate method to approx- 
imate the CDF of a random variable. Knowing k{s) = log [e''''^] - the CGF of 7, the CDF of 7 can be 
approximated by [7, Ch. 1] 



Q{-w) + (j){w) ■ ( 4 - 4 ) if X ^ E^[7] 
w u 



1 fi:'"(0) 



^ 2 6v^[/^"(0)] 



3/2 



(4) 

if X = E^[7] 



with 



w 



sign(s)-y/2(s ■ a; — u = sa/ k"{s), 



where k'{s), k"{s) and k"'{s) are respectively the first, the second and the third derivatives of k{s), 
(f){x) = exp( ), Q{x) = / (j){x)dx, and s depends on x through the saddlepoint equation 



^271 



X 



k'(s) = x; (5) 



solving ^ is the most difficult step of the SPA-based approach. 

In our case, k(s) = Xl/^i i^ii^) once s is found, we can use the relevant expressions for the p-th 
derivative of the CGF of ki{s), shown in TableU to calculate k'{s) = J2l'=i ^'li^) ^^'^ = J^iLi '^'I'i^) 
that are then used to obtain Fy(x) in ©. 

Note that with the SPA approach, 7^ do not have to follow the same fading distribution, i.e., the outage 
probability can be obtained, e.g., when Nakagami-m and Hoyt (or any other combination of non-identical 
distributions) fading signals are combined at the receiver. 

While (HI) is a valid approximation for any value of the argument x, the expression can be further 
simplified noting that we are mostly interested in small values of the outage probability. In particular, 
assuming that x < E^[7] implies that s < (this is because we know that k,"(s) > and k'{s) = 0). So, 



using the well known bound erfc(t) < ^^exp(— t^) |[T4l Appendix 11^ 



we replace Q{—w) by 



ty/n w 

^Or, equivalently, Q(t) < 4>{t)/t. This bound holds for t > and becomes increasingly tight with increasing t. 
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Table II 

Solution of the saddlepoint equation obtained combining L signals with the same distribution. 



and obtain the following approximation of the CDF 

FJx) ^ FJx) = --V^ = , = ) ' (6) 

^ ^ u \sy2Tin"{s) \sy2TXK"{s) 



We thus recover SPA forms similar to those used, for example, in [[14 



A. Identically Distributed Fading 

In the case of mutli-antenna receivers or ARQ transmissions, the channel between the transmitter and 
the receiver is essentially the same for each diversity "branch". Thus, it is reasonable to assume that 7^ 
(/ = 1, . . . , L) are identically distributed. Then, ki{s) = hii(s) and the saddlepoint equation (|5]) reduces to 

L-k[{s)=x. (7) 

In the cases of the fading distributions we consider, from Table |I] we can easily see, that the solution of 
d?]) is obtained solving the quadratic equation. These solutions are shows in Table HIl 

For Nakagami-m fading, the solution is particularly simple and allows us to write the resulting approx- 
imation of the outage probability as 

FJx) = ^ ' ^> exp(-mL ■ ixlL^-l)), (8) 

while the exact solution is known in this case and given by 

F,(x) = r(^mL,^^, (9) 

where r(s,x) = — — / t'^~^e~*dt is the normalized lower incomplete gamma function and T(s) = 

r(s) Jo 

r(s, 00). 



B. Nakagami-m Fading and L = 2 

In the context of opportunistic relay-based communications, an equivalent three-terminal setup (source, 
selected relay, and destination) is quite common [|4j. In such a case, 71 (SNR between the source and 
the destination) and 72 (SNR between the relay and the destination) may have different average SNRs 
7i 7^ 72 when the distances between nodes are different). In this interesting scenario, we can also 
provide a closed-form solution when 71,72 follow Nakagami-m distributions. Equation (|5]) yields 

"^^^^ + =x (10) 



mi — S'ji m2 — S72 

which has the following solution 



i = *_./(]5_™/_!M£^ (U) 

2x y 2x x^ 
with 7 = 7i + 72, 7 = 7i ■ 72, m = mi + m2, m = mi ■ m2, and /3 = (mi/71 + ^2/12) l'^- 



C. General Case 

In the general case, we are not able to solve (|5]) in closed- form. However, the solution of k'{s) = x may 
be found in a few recursive steps which implement the Newton method for solving non-linear equations 

Sk = Sk-i H — — , for /c = 1, . . . , K^^ (12) 

[Sk-l) 

where sq is the initial solution, and the number of recursive steps, -fCmax, controls the accuracy of the 
solution. For the numerical examples we will discuss later, a satisfactory accuracy of the approximation 
is obtained with a relatively small i^max < 5. 

While the requirement of finding s recursively is a drawback of the SPA-based method compared 
to closed-form solutions shown before, we note that the existing altematives often imply a similar 
computational effort. For example, finding the exact form of the CDF F^(x) for Nakagami-m variables 
calls for finding infinite series of coefficients defined recursively [[T]|, a numerical integration implies a 
summation over the integration nodes [|6l, and approximations shown in Q require solving non-linear 
equations. Thus, from the complexity point of view, SPA is a simple approach to calculate the outage. 



IV. Numerical Examples 

In this section, we present a few numerical results to assess the accuracy of the proposed SPA-based 
evaluation of the outage probability. 
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In Fig. [H we show the results when 7 is the sum of L identically distributed 7^, each with 7; = 5dB. 
Three cases are presented: a) Nakagami-m fading signals with nii = 2, b) Rice fading with Ki = 2, c) 
Hoyt fading with qi = 0.5. 

In Fig. [2I we show the results for Nakagami-m distributed signals with L = 2, 71 = 5, 8, 12dB and 
72 = 27]^. We consider two cases for the shape parameters: a) mi = 1 and m2 = 2, and b) mi = 0.5 
and m2 = 2.5. Case (a) can also be solved using the closed-form expressions from [|4l, while case (b) 
requires expressions given in [1] with a truncation of the infinite series. Our method solves both cases 
(a) and (b) in closed-form. 

Finally, we show in Fig. [3] the outage results after combining L = 2 differently fading signals corre- 
sponding to a Nakagami-m fading with parameters m = 1.5, 7i = 5,8, 12dB and a Rice fading with 
K = 5 and 72 = 271. The saddlepoint solution was solved using (fT2l) with -ft'max = 5. 

These examples illustrate well how exact the SFA method is, yielding results that are practically identical 
to the exact form of the CDF. The simplification ^ also provides very accurate results for all outage 
values below 10^^ which, in most cases will be the region of interest when evaluating the performance 
of practical diversity combining schemes. 
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V. Conclusions 

In this work we proposed to use a saddlepoint approximation to evaluate the outage probability 
at receivers combining signals from arbitrarily drawn distributions. We have shown that knowing the 
corresponding moment generating functions, we are able to accurately estimate the outage for an arbitrary 
number of combined signals. The solution requires, in general, solving a scalar non-linear equation whose 
solution can even be found, in particular cases, in closed-form. 
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Figure 3. Exact CDF F^,{x) is compared with the SPA F^{x) and the simpHfied SPA F-y(x) for L = 2 when one variable corresponds to 
Nakagami-m fading with m = 1.5 and another - to the Rice fading with _fL = 5; 72 = 27^. 



